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^r ' We discuss a supersymmetric inverse seesaw model in which lepton flavour violating decays can 



be enhanced either by flavour violating slepton contributions or by the non-unitarity of the charged 
current mixing matrix. As an example we calculate Br{fi -^ cy) taking into account both heavy 
lepton exchange as well as supersymmetric diagrams in a minimal supergravity framework. We 
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expectations, with or without supersymmetry. 

PACS numbers: ll.30.-j, ll.30.Hv, 26.65. +t, 13.15.+g, 14.60.Pq, 95.55.Vj 
Keywords: supersymmetry; neutrino mass and mixing 



* Electronic address: ldepp isch@physik.uni-wuerzburg.de| 
^Electronic address: .valleQific.uv.esI 



1 



I. INTRODUCTION 



;he framework of the standard non-supersymmetric seesaw 
p, y, III, |S{ rare lepton flavour violating decays like // — > 67 are 



A heroic effort dating back to over four decades has finally led to the discovery of neutrino 
oscillations, through a combination of solar, atmospheric, reactor and accelerator neutrino 

npfin 

experiments [H, 12; S 1^ • This has now firmly established the incompleteness of the Standard 
Model of electroweak interactions, suggesting that lepton flavour violation (LFV) may also 
take place in the charged lepton sector. 

It is well known that within t 
models of neutrino masses 
indeed extremely rare, far beyond the sensitivity of any foreseeable experiment. In these 
models the effective light neutrino mass rrii, is inversely proportional to the scale of the lepton 
number violating SU(2) ® U(l) singlet Majorana mass Mr, mu oc M^"*^, while typical LFV 
rates vanish with higher powers of M/j. 

There are two popular mechanisms where lepton flavour violation can be greatly en- 
hanced. First, by supersymmetrizing the model, in which case flavour violation in the 
neutrino Yukawa couplings is automatically transmitted to charged leptons through charged 



slepton and sneutrino loops 
rates. 



y,y 



ll|, giving a sizeable enhancement to the LFV decay 



The other possibility is to consider variants [ij] of the seesaw scheme characterized by 
a small effective lepton number violating Majorana mass term /i so that m,, oc fi. Its 
smallness may be technically natural, since as rrij, — ;► a larger symmetry is achieved |l3|, 
namely lepton number is restored and neutrinos become massless. Due to the fact that in 
this case nii, ^ as fi —>■ this scheme may be called "inverse seesaw". The resulting 
LFV rates do not depend at all on the magnitude of the lepton number violating scale fi, 
which can lie below the weak scale. Moreover, LFV processes may take place even in the 
limit where lepton number is strictly conserved jlj]. Similarly, CP violation can arise in 
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16|. The model thus serves to 



the limit where the light neutrinos are strictly massless 

elucidate that, from a basic point of view, neutrino masses do not play a fundamental role 

in generating flavour violating processes, which in this case are mediated by the exchange 



of SU(2) U(lL singlet leptons, due to the structure of the electroweak charged and 
neutral currents [7j . However, in contrast to the standard seesaw, in inverse seesaw schemes 
the SU(2) ® U(l) singlet leptons need not be super-heavy, leading to highly enhanced 
LFV rates irrespective of the massiveness of neutrinos and irrespective of the existence of 



supersymmetric states 
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This paper is organized as follows. In Sec. Owe introduce the main features of the model, 
while in Sec. Illll we give the renormalization group evolution both of the neutrino sector as 
well as the slepton sector. In Sec. |lV|we discuss the relative importance of the two types of 
contribution to l^ ^ ^Jl decays and in Sec. |V|we give our numerical results and summarize 
the findings in Sec. IVII 



II. INVERSE SEESAW MECHANISM 



^^ii^i^ '-'«) 



(1) 



The particle content of left-handed leptons in the model extends minimally that of the 
Standard Model, by the sequential addition of a pair of two-component SU(2) ® U(l) singlet 
fermions, as follows 

ei 

with i a generation index running over 1,2,3. In addition to the right-handed neutrinos 
characteristic of the standard seesaw model, the inverse seesaw scheme contains an equal 
number of gauge singlet neutrinos Si. In the original formulation of the model, these were 
superstring inspired E{6) singlets, in contrast to the right-handed neutrinos, members of 
the spinorial representation. A similar construction at the SU{2)l ® SU{2)r ® ^(1) level 
was considered in Ref. jl3]- 

The model is characterized by the following symmetric 9x9 mass matrix M. in the iy,iy'^, S 
basis. 



' m^ ^ 



M 



niD M^ 
\ M /i / 

where mo and M are arbitrary 3x3 complex matrices in flavour space, whereas ^ is complex 



(2) 



symmetric. The matrix Ai can be diagonalized by a unitary mixing matrix U^, 

U^M U, = diag(m, , M4 , . . . , Mg) , (3) 

yielding 9 mass eigenstates n^, three of them corresponding to the observed hght neutrinos 

with masses nii, plus the three pairs of two-component leptons {u^, Si) combining to form 

three heavy quasi-Dirac leptons 20|. 

The mass eigenstates are then related to the light neutrino flavour states Ui via the unitary 

matrix 1/^ 

9 

l^i = ^{Uu)iana- (4) 

nn 



which has been studied in earlier papers [IJ, |l5|, |l6|. Assuming ttib, fi <^ M the diagonal- 
ization results in an effective Majorana mass matrix for the light neutrinos hi 

my = m^M'^'^fiM'^mD, (5) 

which may be estimated as 

(Jn^^ = ( "^^ y f-t^) ( ^ ] "' (6) 

VO.leV/ VlOOGeV/ VlkeV/ VlO^GeVy ' ^ ^ 

which vanish in the limit ^ ^ where lepton number conservation is restored. In models 

where lepton number is spontaneously broken by a vacuum expectation value (a) ji21i] one has 

/x = A (a). For typical Yukawas A ~ 10"'^ one sees that /i = 1 keV corresponds to a low scale 

of L violation, (a) ~ 1 MeV. Although such a low scale is not protected by supersymmetry, 

this is not needed, as it is already protected by gauge symmetry. Indeed, being a gauge 

singlet, the smallness of (a) is not destabilized by gauge loops and is technically natural ^. 

In such an "inverse seesaw" scheme the three pairs of singlet neutrinos have masses of 
the order of M and their admixture in the light neutrinos is suppressed as ^. 

In contrast, in the standard seesaw mechanism where the gauge singlet neutrinos Si are 
absent one would have 

I , mD<^Mii^my = m\)Mj^^mD ■ (7) 

niD Mr 



^ For very low values this might lead to interesting signatures in neutrinoless double beta decays [2 



Note that although M and M/j both take the role of a large mass scale suppressing the 
light neutrino masses, their physical meaning is quite different, the former being a Dirac 
mass (AL = 0) and the latter a Majorana mass (AL = 2). In contrast to the mass Mr of 
the right-handed neutrinos of the standard seesaw model, the mass M of our heavy leptons 
can be much smaller, since the suppression in Eq. (jUj) is quadratic and since we have the 
small independent parameter /x characterizing the violation of lepton number. As a result 
the value of M may be as low as the weak scale ^. 

Without loss of generahty one can assume /x to be diagonal, 

fi = diag fii, (8) 

and using the diagonalizing matrix U of the effective light neutrino mass matrix ttt-j,, 

U'^rriiyU = diag rrij, (9) 

equation (jSJ can be written as 



1 = diag ym^ ■ U nijjM ■ diag y/JTi ■ diag ^/Jm ■ M nioU ■ diag y m^ . (10) 

In the basis where the charged lepton Yukawa couplings are diagonal the lepton mixing 
matrix is simply the rectangular matrix formed b y th e first three rows of Uy ^J . 

In analogy to the standard seesaw mechanism \2(y\ it is thus possible to define a complex 
orthogonal matrix 



R = diag ^/J^i ■ M ^moU ■ diag ym^ ^ (11) 

with 6 real parameters. Using R, the neutrino Yukawa coupling matrix Y^, = ^ ^^^ „ rriD can 

be expressed as 

1 



n = —5M ■ diag \Uii^ ■ R ■ diag v^ ■ U\ (12) 

To simplify our discussion we make the assumption that the eigenvalues of both M and 
/i are degenerate and that R is real. This allows us to easily compare our results with 



those obtained previously in Ref. 



271 



28l | for the case of the standard seesaw mechanism. 



^ If light enough, these neutral leptons would be singly-produced in Z^ decays [23|, |2J|, a possibility now 
ruled out by LEP 0. 



The combination Y^LYy which is responsible for flavour non-diagonal slepton mass terms in 



Eqs. (|;^2I;M|) is given by 

{Y}LY^\, = -T^— 1^ %^(^ ■ diag m, ■ U%. (13) 

v^ sm p /U M 

This should be compared with the elements of the light neutrino mass matrix 



[m 



UjlJ 



{U ■ diag m, ■ U^),^ (14) 



The physical consequence of the simplifications we use is that the pattern of LFV transmitted 
to the sleptons is closely correlated to that of the light neutrino sector. The only difference 
being the roles of CP violating phases in (|13|) and (J14j) that may be present in U. In the 
case of CP conservation which we will consider, the correlation is exact. 

In the standard supersymmetric seesaw mechanism the flavour non-diagonal slepton mass 
terms would be (in an analogous approximation, i.e. Mi = Mr and real R matrix) 

(FjLy.),, = ^^^M^ln ^^{U ■ diag m, ■ U%, (15) 

v^ sm p Mfj 

as compared to Eq. UHl 

III. RENORMALIZATION GROUP EVOLUTION 
A. The neutrino sector 

In what follows it will be sufficient for us to confine ourselves to the simpler case where 
heavy neutrino masses are degenerate. Below the scale of M, the one-loop renormalization 
group equation (RGE) for the effective neutrino mass matrix in the MSSM is given by 



d 

dt 16 



b (("^^'' " l^' + Tr{6YiYu)] m, + ((i;^l^)m, + m,{Y!Yif)\ , (16) 



with the U(l) and SU(2) gauge couplings gi and (72, and the Yukawa coupling matrices Yj/ 
and Yi for the charge |-quarks and charged leptons, respectively ^. This RGE is linear in 



The corresponding evolution equations for gi,2, Yjj and Y/ can be found in 



H 



m^ and can thus be solved analytically 



as 



m^it) = mm^m, t = In (^^ . (17) 

Since the evolution is dominated by the gauge and third generation Yukawa couplings one 
obtains, to a good approximation: 

J(t) = Vi-diag(l,l,/,) (18) 

with 

W = exp (^ j\-3gl - ^gl)dt'^ , (19) 



1 



^*w = e^p(T^y^ mi'dt'], (20) 

1 '" 



Irit) =^^P[y^2 1 \yr\'dt'j. (21) 

Above the scale M, which can be substantially lower than the corresponding heavy lepton 
scale in standard seesaw schemes, the evolution of the neutrino Yukawa matrix Y^ is governed 

by y 

j/" = T^^^ (("^^' " ¥' + '^'^(^^uYu + YX)^ l + Y,^Yi + 3YX^ ■ (22) 

B. The slepton sector 

Having evolved the neutrino Yukawa couplings to the unification scale Mquti one can now 
run the slepton mass matrix from Mqut to the electroweak scale assuming the mSUGRA 
universality conditions on the soft SUSY breaking terms m| (left-handed slepton doublets), 
m\ (charged right-handed slepton singlets), m^ (right-handed sneutrino singlets), Ae (tri- 
linear couplings analogous to Ye) and A^ (trilinear couplings analogous to Yi,) at Mqut'- 

ml = ml = m% = mll, A^ = A^Y^, A, = AoY„ (23) 

where ttiq is the common scalar mass and Aq the common trilinear coupling. For definiteness 
we adopt in our following analysis, the mSUGRA benchmark scenario SPSla proposed in 



3l|, described by 

mo = 100 GeV, mi/2 = 250 GeV, A = -100 GeV, tan /3 = 10, /i > 0, (24) 

where // is here the SUSY Higgs-mixing parameter. In general, the charged slepton (mass)^ 
matrix has the form: 

.^J< ^J \ , (25) 

where m? , m? and m? are 3x3 matrices, m? and m? being hermitian. The matrix 

IL IR ILR IL IR 

elements are given by 

(m?J„5 = {ml)ab + 5ab Iml + m|cos(2/?) f-- + sin'^^y j j (26) 

{m^Jab = {rn%)ab + Sab{m^^ - m% cos (2/?) sin^ Ow) (27) 

{nn?r )ab = {A^) abV cos (3 - dabmi^f^ tan f3. (28) 



Applying the mSUGRA conditions (j23|) at Mqut and performing the evolution to Mz the 
SUSY breaking terms can be expressed as: 

m\ = mll + ((5m|)MssM + Sml (29) 

mji = mil + {6m%)MssM + 5m\ (30) 

Ae = AoYe + Mmssm + SA, (31) 

with ((5m| ^)mssm and ((5v4)mssm denoting the usual flavour diagonal MSSM renormalization 
group corrections 30|. In addition, the presence of right-handed neutrinos radiatively in- 
duces flavour off-diagonal terms denoted by Smj^ ^ and 6 A. In the leading-log approximation 
these terms are given as 3j| 

^< = -^{3rnl + Al){Y^LY,) (32) 

6m% = (33) 

SA = -^ye ■ (Y^LY,) (34) 



with 



L = diag(l„^). (35) 



Finally, the physical charged slepton masses are then found by diagonalizing ()25|) using the 
6x6 unitary matrix Up 

^m} Ui = diag(m? , • • • , m| ) . (36) 

Correspondingly, the slepton mass eigenstates are expressed in terms of the gauge eigenstates 
by 

L = {Uf)ialLi + {U^){i+3)jRi, a = 1, ..., 6; i = e,fi,T. (37) 

Similarly to (j26|) . the 6x6 (mass)^ matrix of the SUSY partners of the left- and right-handed 
neutrinos is given by 

-nffi-'l. (38) 

\ ^LR ^R / 

with 

(mljab = (mDab + - 6 abml cos{2P) + {m^)li, (39) 

(mljab = im%U + ^2, (40) 

("^^ifl)afe = V COS (3{A^)ab - 1^ cot P{mD)ab- (41) 

The sneutrino mass matrix is diagonalized by a unitary 6x6 matrix [/;>, 

UlmjU^ = diag(m?^, ...,m?J, (42) 

which, like f/j,, can contain sizeable mixings between SU(2) U(l) isodoublet and isosinglet 
sneutrinos. In contrast to the standard seesaw where the right-handed neutrinos (and thus 
sneutrinos) are extremely heavy and can be safely neglected in low energy processes, this 
is not the case in our "inverse seesaw" model. Indeed, the presence of the gauge singlet 
superfields Si is crucial in this model in causing a big enhancement in the LFV rates. Note 
however that the scalars present in Eq. (jHHj) can be neglected in calculating the LFV decay 
rates. 

IV. LEPTON FLAVOUR VIOLATING DECAYS: Ir -^ Ij^ 

The effective Lagrangian for l^ -^ Ij'-f may be written generically as 

Ceff = Iha^pF-' {AlP, + AlPn) k, (43) 



where F°^ is the electromagnetic field strength tensor, aap = | [7a, 7/3] and Prl = |(1 ±75) 
are the helicity projection operators. The coefficients A^l ^ are determined by the relevant 
gauge theory Fe ynm an diagrams. The decay rate for l^ -^ Ij'^ that follows from PHj) can 
be expressed as J22| 

r(c-^77)=f^f.(i^^r+i4r)- (44) 



A. Heavy lepton contribution 



The contribution to the decay /^ —>■ L 7 arising from the admixture o: 



Q 



the heavy neutrinos 



in the left-handed charged current SU(2) ® U(l) weak interaction [7| exists both for the 
seesaw scheme as well as for the inverse seesaw model. One finds that the branching ratio 

m 



IS given as 






where 

r2 



G'.,=2^(f/.):,(f/.),fcG,(^^J, G,{x) = -^-_-^--^_-^lnx, (46) 



k 

Fj is the total decay rate of lepton i and U^ is the matrix describing the diagonalization of 
the neutrino mass matrix in the seesaw scheme under consideration. As noted in [lj| within 
the framework of the inverse SU(2) ® U(l) seesaw mechanism the decay /~ -^ /~7 occurs 
in the limit of lepton number conservation where the light neutrinos become massless. The 
rate is correspondingly enhanced with respect to that of the simplest seesaw scheme, as it 
is not suppressed by the smallness of neutrino masses. 

B. Supersymmetric contribution 

The coefficients v4^^ of the effective Lagrangian for l^ —>■ Ij'-f in Eq. ()43p in the MSSM 



have been given in Ref. [35[. They are determined by the photon penguin diagrams shown 
in Fig. n with charginos/sneutrinos or neutralinos/charged sleptons circulating in the loop. 
The superscript c (n) refers to the chargino (neutralino) diagram of Fig. ^ while the ffavour 
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Figure 1: Supersymmetric diagrams for l^ — > L- 7 

indices are omitted. Because m/. ^ mi. and m? is diagonal (see ()27I30I33|) ). one has 

. The dominant amphtudes in ()44p are approximately given by 

2 



A% 



1 



gimi. 



327r2 y/2mw cos /3 



EE 



m 



Xa 



X- 



-3 + Kk 



Tfl- 
a=l k=l ^k 



21nr^fc 



{On)al{OL)a2{U;),k{U,),k 



^"aU? 



ak) 



(47) 



A'r 



1 



6 



-gltunew^^ 



Act 



327r2 



2 (0Ar)al((0Ar)a2+ (0Ar)aitan6'H/) 
a=l k=l Ik 



' ak) 



with 



ak 



m- 



m 



^k 



ak 



m^o 
mr 



(48) 



(49) 



the chargino diagonalization matrices Ol,Or and the neutralino diagonalization matrix 
Oat. The mass eigenvalues of the charginos and neutralinos are denoted by m-- and m^o, 
respectively. It is worth noting that the sum in (|47|) runs over both left- and right-handed 
sneutrinos. The right-handed sneutrinos will yield a sizeable contribution for small M 
whereas the left-handed ones and the charged sleptons become significant for larger M. 



The contribution due to flavour non-diagonal slepton mass terms is roughly given by 

2 



Br{li -^ Ij-f) ^ a^ tan^ /? 



T^i \ 'mi. 

''I \ '•I 



[5m 



LJij 



m? 



(50) 



m / Ti 
where m is the approximate mass scale of SUSY particles in the loops. 

The numerical calculations discussed later are performed with the full expressions for 



A^" and A^", which can be found in 



33I and lag 
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Figure 2: Contours of Br(fj, — > 67) in the (M, -^)-plane (logarithmic scales) for hierarchical light 
neutrinos with mi = eV (left panel) and for degenerate light neutrinos with m-i = 0.3 eV (right 
panel). The dark (light) area is excluded for Br{fi — > 67) < 10^^^ (10^^^). The blue contours in 
the lower left depict the contribution from neutral heavy leptons only. The diagonal lines show 
contours of constant /^ = 1, 10^'^, 10^^, 10~^ GeV (top to bottom). The vertical lines are contours 
of Br[^ — > 67) in the standard SUSY seesaw (see text). 

V. NUMERICAL RESULTS 

Low energy neutrino experiments |l|, |2; ld| now provide substantial information on the 
light neutrino masses and lepton mixing matrix. This information has to be evolved to the 
unification scale in order to calculate the slepton mass corrections. Our numerical calculation 

n 

is performed as follows. We fix the light neutrino sector by using the latest global fit |^ 
for the neutrino oscillation parameters, neglecting possible CP phases to which current data 
are insensitive. The Yukawa coupling Y^ is then calculated via Eq. p2|) . The result for 
Br{^ -^ 67) in the case of hierarchical neutrino masses, mi = eV, is shown in Fig. |21 (left 
panel) as a contour plot in the (M, -p)-plane, with /i < O.IM on the whole plane ^. The dark 
(red) area is excluded by the current experimental limit Br{fi — > 67) < 1.2 ■ 10^^^ while the 



The diagonal lines depict contours of constant /i in the inverse seesaw case, for easier reading. 
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light (brown) band shows the additional sensitivity aimed at the PSI experiment Br{fi — >■ 
67) < 1.5 ■ 10^^^. For large M, the RG&induced flavour violating slepton contributions are 
clearly dominant due to the i -jjj j In^ ^^^^"^ dependence from Eqs. (IT^ and (fSUj) . It is only 
below M ^ 10^ GeV where the admixture of the heavy singlet neutrino (and sneutrino) 
states becomes significant and dominant. In order to guide the eye we have also included 
the corresponding contours for the standard SUSY seesaw // ^ 67 decay branching ratio 
at 10"^^ and 10~^^ (dashed), that would arise from taking Mr = M, as indicated by the 
vertical lines on the right side. These are excatly vertical, as the seesaw does not contain the 
parameter fi. We see that for given Mr = M value the rate can be enhanced with respect 
to seesaw model expectations, with or without supersymmetry. 

The curved (blue) contours peaked at slightly above M=100 GeV correspond to the non- 



M 



liU- 



The 



supersymmetric isosinglet neutral heavy lepton contribution, considered in 
difference is that now we take into account the neutrino masses indicated by neutrino os- 
cillation data. This actually has no impact, as the violation of flavour in this case arises 
mainly from the isosinglet neutral heavy lepton contribution, which in this model is essen- 
tially unrelated to the light neutrino masses, due to the freedom in choosing the value of fi. 
On the other hand, in contrast to the above works we explicitly correlate the LFV decay 
rate to the neutrino oscillation mixing angles. Finally, the rise on the very far left of the 
plots corresponds the the contribution in the supersymmetric loops arising from SU(2) ® 
U(l) singlet scalar neutrinos. 

The right panel in Fig. |21 shows the analogous plot for quasi-degenerate neutrino masses, 
rrti = 0.3 eV, where one can see that Br{fi — * 67) is suppressed by roughly one order of 
magnitude as compared to the hierarchical neutrino case. This is because the flavour non- 
diagonal elements of the light neutrino mixing matrix (which are ultimately the source of 
LFV, both RGE- and non-universality-induced) are suppressed as — — - for large mi 
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VI. CONCLUSIONS 

We have discussed a supersymmetric inverse seesaw model in which lepton flavour vio- 
lating decays can be enhanced either by flavour violating slepton contributions or due to 
heavy lepton exchange. In contrast to the standard seesaw scheme, the SU(2) ® U(l) isos- 
inglet heavy leptons present in this model can be relatively light and contribute significantly 
to LFV processes, irrespective of the magnitude of neutrino masses and irrespective of the 
supersymmetric contributions. We have considered both types of contributions focusing on 
their differences with respect to the standard supersymmetric seesaw scheme. As an exam- 
ple we have calculated Br{fi -^ 67) taking into account both neutral heavy lepton as well 
as supersymmetric diagrams in a minimal supergravity framework. Clearly, additional LFV 
processes such as fj, ^ 3e, fi ^ e conversion in nuclei, as well as LFV r decays can be 
considered generalizing to the supersymmetric case the analysis presented in 17|, |3J| . 
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